{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}C:\07102020\research, journal papers\nationalism and ethnic terror\do and log-files\replication materials\do and log-files\national_terror, log-file.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}12 Mar 2021, 18:51:42

{com}. do "C:\07102020\research, journal papers\nationalism and ethnic terror\do and log-files\replication materials\do and log-files\national_terror, do-file.do"
{txt}
{com}. * March 12, 2021
. * national_terror, do-file
. 
. set matsize 500
{txt}
{com}. set more off
{txt}
{com}. 
. use "C:\07102020\research, journal papers\nationalism and ethnic terror\do and log-files\replication materials\data\national_terror, data.dta", clear
{txt}
{com}.          
. 
. ***** Table 1 *****
. sum attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count ///                                                                                                           
>         if isrelevant==1 

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 6}attack {c |}{res}     22,754    1.021403    9.467018          0        317
{txt}v2exl_legi~d {c |}{res}     22,754    .5018768    .3095542          0          1
{txt}{space 4}xpolity2 {c |}{res}     22,754    2.046629    3.835802         -6          7
{txt}status_aut~y {c |}{res}     22,754    .1630922    .3694579          0          1
{txt}status_pow~s {c |}{res}     22,754    .2895755    .4535753          0          1
{txt}{hline 13}{c +}{hline 57}
status_dis~m {c |}{res}     22,754     .136899    .3437483          0          1
{txt}status_sep~t {c |}{res}     22,754    .0172717    .1302848          0          1
{txt}Weighted_C~g {c |}{res}     22,754    4.498575    3.957805          0    10.8638
{txt}{space 3}cw_binary {c |}{res}     22,754    .2595148    .4383781          0          1
{txt}log_rgdppc~g {c |}{res}     22,754    8.091533    1.231858   4.764052   11.48176
{txt}{hline 13}{c +}{hline 57}
muspct_lag~d {c |}{res}     22,754    24.23301    33.37175          0   99.73901
{txt}{space 3}groupsize {c |}{res}     22,754    .1754934      .25724      .0001        .98
{txt}excl_group~t {c |}{res}     22,754    10.48958    16.62133          0         55
{txt}
{com}. 
. 
. ***** Table 2 *****
. * Model 1
. menbreg attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count ///                                                                                                           
>         if isrelevant==1||cowcode: ||cowgroupid: ,cov(uns)
{res}{txt}
Fitting fixed-effects model:

Iteration 0:{space 3}log likelihood = {res:-26924.106}  
Iteration 1:{space 3}log likelihood = {res:-14798.912}  
Iteration 2:{space 3}log likelihood = {res:-9704.0422}  (not concave)
Iteration 3:{space 3}log likelihood = {res:-9430.1177}  (not concave)
Iteration 4:{space 3}log likelihood = {res:-9090.0742}  (not concave)
Iteration 5:{space 3}log likelihood = {res:-8906.7726}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-8727.8768}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-8639.3288}  (not concave)
Iteration 8:{space 3}log likelihood = {res:-8581.2962}  (not concave)
Iteration 9:{space 3}log likelihood = {res:-8526.7546}  
Iteration 10:{space 2}log likelihood = {res: -8477.845}  
Iteration 11:{space 2}log likelihood = {res:-8166.0768}  
Iteration 12:{space 2}log likelihood = {res:-8149.4073}  
Iteration 13:{space 2}log likelihood = {res:-8149.2749}  
Iteration 14:{space 2}log likelihood = {res:-8149.2748}  

Refining starting values:

Grid node 0:{space 3}log likelihood = {res:-7722.2707}

Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-7722.2707}  (not concave)
Iteration 1:{space 3}log likelihood = {res:-6744.3902}  
Iteration 2:{space 3}log likelihood = {res:-6562.1608}  
Iteration 3:{space 3}log likelihood = {res:-6480.4007}  
Iteration 4:{space 3}log likelihood = {res:-6454.2568}  
Iteration 5:{space 3}log likelihood = {res:-6448.4219}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-6447.9885}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-6448.3029}  (not concave)
Iteration 8:{space 3}log likelihood = {res:-6448.3952}  
Iteration 9:{space 3}log likelihood = {res:-6448.4031}  
Iteration 10:{space 2}log likelihood = {res:-6448.2592}  
Iteration 11:{space 2}log likelihood = {res:-6448.1206}  
Iteration 12:{space 2}log likelihood = {res:-6448.1247}  
Iteration 13:{space 2}log likelihood = {res:-6448.1781}  
Iteration 14:{space 2}log likelihood = {res:-6448.2352}  
Iteration 15:{space 2}log likelihood = {res:-6448.2823}  
Iteration 16:{space 2}log likelihood = {res:-6448.3181}  
Iteration 17:{space 2}log likelihood = {res:-6448.3446}  
Iteration 18:{space 2}log likelihood = {res: -6448.364}  
Iteration 19:{space 2}log likelihood = {res:-6448.3782}  
Iteration 20:{space 2}log likelihood = {res:-6448.3886}  
Iteration 21:{space 2}log likelihood = {res:-6448.3963}  
Iteration 22:{space 2}log likelihood = {res:-6448.4018}  
Iteration 23:{space 2}log likelihood = {res:-6448.4059}  
Iteration 24:{space 2}log likelihood = {res:-6448.4088}  
Iteration 25:{space 2}log likelihood = {res: -6448.411}  
Iteration 26:{space 2}log likelihood = {res:-6448.4126}  
Iteration 27:{space 2}log likelihood = {res:-6448.4137}  
Iteration 28:{space 2}log likelihood = {res:-6448.4146}  
Iteration 29:{space 2}log likelihood = {res:-6448.4152}  
{res}
{txt}Mixed-effects nbinomial regression{col 49}{txt}Number of obs{col 67}={res}{col 69}    22,754
{txt}Overdispersion: {col 28}{res}mean

{txt}{hline 16}{c TT}{hline 44}
{col 17}{txt}{c |}{col 23}No. of{col 36}Observations per Group
{col 2}{txt}Group Variable{col 17}{c |}{col 23}Groups{col 33}Minimum{col 44}Average{col 55}Maximum
{txt}{hline 16}{c +}{hline 44}
{col 9}{res}cowcode{col 17}{txt}{c |}{res}{col 21}     132{col 31}        6{col 42}    172.4{col 53}    1,899
{col 6}{res}cowgroupid{col 17}{txt}{c |}{res}{col 21}     757{col 31}        1{col 42}     30.1{col 53}       40
{txt}{hline 16}{c BT}{hline 44}

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}12{txt}){col 67}={res}{col 70}   527.05
{txt}Log likelihood = {res}-6448.4152{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2} 1.008113{col 45}{space 2} .3123017{col 56}{space 1}    3.23{col 65}{space 3}0.001{col 73}{space 4} .3960127{col 86}{space 3} 1.620213
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0461254{col 45}{space 2} .0180871{col 56}{space 1}    2.55{col 65}{space 3}0.011{col 73}{space 4} .0106753{col 86}{space 3} .0815754
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2} .3536592{col 45}{space 2} .3026589{col 56}{space 1}    1.17{col 65}{space 3}0.243{col 73}{space 4}-.2395413{col 86}{space 3} .9468597
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} .0573547{col 45}{space 2} .2223731{col 56}{space 1}    0.26{col 65}{space 3}0.796{col 73}{space 4}-.3784887{col 86}{space 3}  .493198
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2} .9868245{col 45}{space 2} .2370772{col 56}{space 1}    4.16{col 65}{space 3}0.000{col 73}{space 4} .5221618{col 86}{space 3} 1.451487
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} 1.203089{col 45}{space 2} .4336893{col 56}{space 1}    2.77{col 65}{space 3}0.006{col 73}{space 4} .3530738{col 86}{space 3} 2.053104
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2} .1069338{col 45}{space 2} .0141784{col 56}{space 1}    7.54{col 65}{space 3}0.000{col 73}{space 4} .0791447{col 86}{space 3} .1347229
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.465608{col 45}{space 2} .1344052{col 56}{space 1}   18.34{col 65}{space 3}0.000{col 73}{space 4} 2.202179{col 86}{space 3} 2.729038
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2} .7421893{col 45}{space 2} .1341913{col 56}{space 1}    5.53{col 65}{space 3}0.000{col 73}{space 4} .4791791{col 86}{space 3}   1.0052
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2} .0115971{col 45}{space 2} .0056515{col 56}{space 1}    2.05{col 65}{space 3}0.040{col 73}{space 4} .0005204{col 86}{space 3} .0226738
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.5291519{col 45}{space 2} .8143785{col 56}{space 1}   -0.65{col 65}{space 3}0.516{col 73}{space 4}-2.125304{col 86}{space 3} 1.067001
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0446202{col 45}{space 2} .0175221{col 56}{space 1}   -2.55{col 65}{space 3}0.011{col 73}{space 4}-.0789629{col 86}{space 3}-.0102776
{txt}{space 26}_cons {c |}{col 33}{res}{space 2}-17.54124{col 45}{space 2} 1.418181{col 56}{space 1}  -12.37{col 65}{space 3}0.000{col 73}{space 4}-20.32083{col 86}{space 3}-14.76166
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 23}/lnalpha {c |}{col 33}{res}{space 2} 1.608561{col 45}{space 2} .0415146{col 56}{space 1}   38.75{col 65}{space 3}0.000{col 73}{space 4} 1.527193{col 86}{space 3} 1.689928
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode                        {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 8.12e-19{col 45}{space 2} 7.00e-10{col 73}{space 4}        .{col 86}{space 3}        .
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode>cowgroupid             {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 41.17977{col 45}{space 2} 6.871314{col 73}{space 4} 29.69281{col 86}{space 3} 57.11059
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{txt}LR test vs. nbinomial model: {txt}chi2({res}2{txt}) ={res} 3401.72{col 59}{txt}Prob > chi2 ={res}{col 73}0.0000

{txt}{p 0 6 4 79}Note: {help j_mixedlr##|_new:LR test is conservative} and provided only for reference.{p_end}

{com}. * Model 2
. sort cowgroupid year 
{txt}
{com}. xtabond2 attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///       
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count L.attack ///                                                                                                           
>         if isrelevant==1, gmm (v2exl_legitideolcr_0_lagged, lag (1 5)) ///
>         ivstyle(xpolity2   status_autonomy status_powerless status_discrim status_separatist ///       
>     Weighted_Conflict_Index_lag_log cw_binary   ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count L.attack)
{res}{txt}Favoring space over speed. To switch, type or click on {stata "mata: mata set matafavor speed, perm"}.

Dynamic panel-data estimation, one-step system GMM
{hline 78}
Group variable: {res}cowgroupid{txt}{col 49}Number of obs      = {res}    22033
{txt}Time variable : {res}year{txt}{col 49}Number of groups   = {res}      756
{txt}Number of instruments = {res}236{col 49}{txt}Obs per group: min = {res}        1
{txt}Wald chi2({res}13{txt}){col 15}= {res}  6580.60{col 64}{txt}avg = {res}    29.14
{txt}Prob > chi2{col 15}={res}     0.000{col 64}{txt}max = {res}       39
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2} 1.841913{col 45}{space 2} .5655657{col 56}{space 1}    3.26{col 65}{space 3}0.001{col 73}{space 4} .7334248{col 86}{space 3} 2.950401
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0351976{col 45}{space 2} .0134742{col 56}{space 1}    2.61{col 65}{space 3}0.009{col 73}{space 4} .0087886{col 86}{space 3} .0616066
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2} 1.283306{col 45}{space 2} .1677924{col 56}{space 1}    7.65{col 65}{space 3}0.000{col 73}{space 4} .9544391{col 86}{space 3} 1.612173
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} .8389213{col 45}{space 2}  .117725{col 56}{space 1}    7.13{col 65}{space 3}0.000{col 73}{space 4} .6081845{col 86}{space 3} 1.069658
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2} 1.195106{col 45}{space 2} .1420869{col 56}{space 1}    8.41{col 65}{space 3}0.000{col 73}{space 4} .9166206{col 86}{space 3} 1.473591
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2}-.1611778{col 45}{space 2}  .369556{col 56}{space 1}   -0.44{col 65}{space 3}0.663{col 73}{space 4}-.8854943{col 86}{space 3} .5631387
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2} .0652949{col 45}{space 2} .0131301{col 56}{space 1}    4.97{col 65}{space 3}0.000{col 73}{space 4} .0395605{col 86}{space 3} .0910294
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 1.581347{col 45}{space 2} .1247306{col 56}{space 1}   12.68{col 65}{space 3}0.000{col 73}{space 4}  1.33688{col 86}{space 3} 1.825815
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2} .2855795{col 45}{space 2} .0381654{col 56}{space 1}    7.48{col 65}{space 3}0.000{col 73}{space 4} .2107768{col 86}{space 3} .3603823
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2}-.0079017{col 45}{space 2} .0014297{col 56}{space 1}   -5.53{col 65}{space 3}0.000{col 73}{space 4} -.010704{col 86}{space 3}-.0050995
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2} .2566621{col 45}{space 2} .1890627{col 56}{space 1}    1.36{col 65}{space 3}0.175{col 73}{space 4} -.113894{col 86}{space 3} .6272181
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0175138{col 45}{space 2} .0059574{col 56}{space 1}   -2.94{col 65}{space 3}0.003{col 73}{space 4}-.0291902{col 86}{space 3}-.0058374
{txt}{space 31} {c |}
{space 25}attack {c |}
{space 28}L1. {c |}{col 33}{res}{space 2} .3588001{col 45}{space 2} .0051531{col 56}{space 1}   69.63{col 65}{space 3}0.000{col 73}{space 4} .3487003{col 86}{space 3}    .3689
{txt}{space 31} {c |}
{space 26}_cons {c |}{col 33}{res}{space 2}-3.617014{col 45}{space 2} .4805071{col 56}{space 1}   -7.53{col 65}{space 3}0.000{col 73}{space 4}-4.558791{col 86}{space 3}-2.675237
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Instruments for first differences equation
{col 3}Standard
{col 5}D.(xpolity2 status_autonomy status_powerless status_discrim
{col 5}status_separatist Weighted_Conflict_Index_lag_log cw_binary log_rgdppc_lag
{col 5}muspct_lagged groupsize excl_groups_count L.attack)
{col 3}GMM-type (missing=0, separate instruments for each period unless collapsed)
{col 5}L(1/5).v2exl_legitideolcr_0_lagged
Instruments for levels equation
{col 3}Standard
{col 5}xpolity2 status_autonomy status_powerless status_discrim status_separatist
{col 5}Weighted_Conflict_Index_lag_log cw_binary log_rgdppc_lag muspct_lagged
{col 5}groupsize excl_groups_count L.attack
{col 5}_cons
{col 3}GMM-type (missing=0, separate instruments for each period unless collapsed)
{col 5}D.v2exl_legitideolcr_0_lagged
{hline 78}
Arellano-Bond test for AR(1) in first differences:{col 52}z = {res}-82.63{txt}  Pr > z = {res} 0.000
{txt}Arellano-Bond test for AR(2) in first differences:{col 52}z = {res} 22.27{txt}  Pr > z = {res} 0.000
{txt}{hline 78}
Sargan test of overid. restrictions: chi2({res}222{txt}){col 49}={res} 954.41{txt}{col 59}Prob > chi2 = {res} 0.000
{txt}  (Not robust, but not weakened by many instruments.)

Difference-in-Sargan tests of exogeneity of instrument subsets:
  GMM instruments for levels
    Sargan test excluding group:{col 38}chi2({res}183{txt}){col 49}={res} 543.69{txt}{col 59}Prob > chi2 = {res} 0.000
{txt}    Difference (null H = exogenous):{col 38}chi2({res}39{txt}){col 49}={res} 410.72{txt}{col 59}Prob > chi2 = {res} 0.000
{txt}  iv(xpolity2 status_autonomy status_powerless status_discrim status_separatist Weighted_Conflict_Index_lag_log cw_binary log_rgdppc_lag muspct_lagged groupsize excl_groups_count L.attack)
    Sargan test excluding group:{col 38}chi2({res}210{txt}){col 49}={res} 616.53{txt}{col 59}Prob > chi2 = {res} 0.000
{txt}    Difference (null H = exogenous):{col 38}chi2({res}12{txt}){col 49}={res} 337.88{txt}{col 59}Prob > chi2 = {res} 0.000

{txt}
{com}. 
.         
. ***** Figure 1 *****
. menbreg attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count ///                                                                                                           
>         if isrelevant==1||cowcode: ||cowgroupid: ,cov(uns)
{res}{txt}
Fitting fixed-effects model:

Iteration 0:{space 3}log likelihood = {res:-26924.106}  
Iteration 1:{space 3}log likelihood = {res:-14798.912}  
Iteration 2:{space 3}log likelihood = {res:-9704.0422}  (not concave)
Iteration 3:{space 3}log likelihood = {res:-9430.1177}  (not concave)
Iteration 4:{space 3}log likelihood = {res:-9090.0742}  (not concave)
Iteration 5:{space 3}log likelihood = {res:-8906.7726}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-8727.8768}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-8639.3288}  (not concave)
Iteration 8:{space 3}log likelihood = {res:-8581.2962}  (not concave)
Iteration 9:{space 3}log likelihood = {res:-8526.7546}  
Iteration 10:{space 2}log likelihood = {res: -8477.845}  
Iteration 11:{space 2}log likelihood = {res:-8166.0768}  
Iteration 12:{space 2}log likelihood = {res:-8149.4073}  
Iteration 13:{space 2}log likelihood = {res:-8149.2749}  
Iteration 14:{space 2}log likelihood = {res:-8149.2748}  

Refining starting values:

Grid node 0:{space 3}log likelihood = {res:-7722.2707}

Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-7722.2707}  (not concave)
Iteration 1:{space 3}log likelihood = {res:-6744.3902}  
Iteration 2:{space 3}log likelihood = {res:-6562.1608}  
Iteration 3:{space 3}log likelihood = {res:-6480.4007}  
Iteration 4:{space 3}log likelihood = {res:-6454.2568}  
Iteration 5:{space 3}log likelihood = {res:-6448.4219}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-6447.9885}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-6448.3029}  (not concave)
Iteration 8:{space 3}log likelihood = {res:-6448.3952}  
Iteration 9:{space 3}log likelihood = {res:-6448.4031}  
Iteration 10:{space 2}log likelihood = {res:-6448.2592}  
Iteration 11:{space 2}log likelihood = {res:-6448.1206}  
Iteration 12:{space 2}log likelihood = {res:-6448.1247}  
Iteration 13:{space 2}log likelihood = {res:-6448.1781}  
Iteration 14:{space 2}log likelihood = {res:-6448.2352}  
Iteration 15:{space 2}log likelihood = {res:-6448.2823}  
Iteration 16:{space 2}log likelihood = {res:-6448.3181}  
Iteration 17:{space 2}log likelihood = {res:-6448.3446}  
Iteration 18:{space 2}log likelihood = {res: -6448.364}  
Iteration 19:{space 2}log likelihood = {res:-6448.3782}  
Iteration 20:{space 2}log likelihood = {res:-6448.3886}  
Iteration 21:{space 2}log likelihood = {res:-6448.3963}  
Iteration 22:{space 2}log likelihood = {res:-6448.4018}  
Iteration 23:{space 2}log likelihood = {res:-6448.4059}  
Iteration 24:{space 2}log likelihood = {res:-6448.4088}  
Iteration 25:{space 2}log likelihood = {res: -6448.411}  
Iteration 26:{space 2}log likelihood = {res:-6448.4126}  
Iteration 27:{space 2}log likelihood = {res:-6448.4137}  
Iteration 28:{space 2}log likelihood = {res:-6448.4146}  
Iteration 29:{space 2}log likelihood = {res:-6448.4152}  
{res}
{txt}Mixed-effects nbinomial regression{col 49}{txt}Number of obs{col 67}={res}{col 69}    22,754
{txt}Overdispersion: {col 28}{res}mean

{txt}{hline 16}{c TT}{hline 44}
{col 17}{txt}{c |}{col 23}No. of{col 36}Observations per Group
{col 2}{txt}Group Variable{col 17}{c |}{col 23}Groups{col 33}Minimum{col 44}Average{col 55}Maximum
{txt}{hline 16}{c +}{hline 44}
{col 9}{res}cowcode{col 17}{txt}{c |}{res}{col 21}     132{col 31}        6{col 42}    172.4{col 53}    1,899
{col 6}{res}cowgroupid{col 17}{txt}{c |}{res}{col 21}     757{col 31}        1{col 42}     30.1{col 53}       40
{txt}{hline 16}{c BT}{hline 44}

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}12{txt}){col 67}={res}{col 70}   527.05
{txt}Log likelihood = {res}-6448.4152{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2} 1.008113{col 45}{space 2} .3123017{col 56}{space 1}    3.23{col 65}{space 3}0.001{col 73}{space 4} .3960127{col 86}{space 3} 1.620213
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0461254{col 45}{space 2} .0180871{col 56}{space 1}    2.55{col 65}{space 3}0.011{col 73}{space 4} .0106753{col 86}{space 3} .0815754
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2} .3536592{col 45}{space 2} .3026589{col 56}{space 1}    1.17{col 65}{space 3}0.243{col 73}{space 4}-.2395413{col 86}{space 3} .9468597
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} .0573547{col 45}{space 2} .2223731{col 56}{space 1}    0.26{col 65}{space 3}0.796{col 73}{space 4}-.3784887{col 86}{space 3}  .493198
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2} .9868245{col 45}{space 2} .2370772{col 56}{space 1}    4.16{col 65}{space 3}0.000{col 73}{space 4} .5221618{col 86}{space 3} 1.451487
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} 1.203089{col 45}{space 2} .4336893{col 56}{space 1}    2.77{col 65}{space 3}0.006{col 73}{space 4} .3530738{col 86}{space 3} 2.053104
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2} .1069338{col 45}{space 2} .0141784{col 56}{space 1}    7.54{col 65}{space 3}0.000{col 73}{space 4} .0791447{col 86}{space 3} .1347229
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.465608{col 45}{space 2} .1344052{col 56}{space 1}   18.34{col 65}{space 3}0.000{col 73}{space 4} 2.202179{col 86}{space 3} 2.729038
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2} .7421893{col 45}{space 2} .1341913{col 56}{space 1}    5.53{col 65}{space 3}0.000{col 73}{space 4} .4791791{col 86}{space 3}   1.0052
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2} .0115971{col 45}{space 2} .0056515{col 56}{space 1}    2.05{col 65}{space 3}0.040{col 73}{space 4} .0005204{col 86}{space 3} .0226738
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.5291519{col 45}{space 2} .8143785{col 56}{space 1}   -0.65{col 65}{space 3}0.516{col 73}{space 4}-2.125304{col 86}{space 3} 1.067001
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0446202{col 45}{space 2} .0175221{col 56}{space 1}   -2.55{col 65}{space 3}0.011{col 73}{space 4}-.0789629{col 86}{space 3}-.0102776
{txt}{space 26}_cons {c |}{col 33}{res}{space 2}-17.54124{col 45}{space 2} 1.418181{col 56}{space 1}  -12.37{col 65}{space 3}0.000{col 73}{space 4}-20.32083{col 86}{space 3}-14.76166
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 23}/lnalpha {c |}{col 33}{res}{space 2} 1.608561{col 45}{space 2} .0415146{col 56}{space 1}   38.75{col 65}{space 3}0.000{col 73}{space 4} 1.527193{col 86}{space 3} 1.689928
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode                        {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 8.12e-19{col 45}{space 2} 7.00e-10{col 73}{space 4}        .{col 86}{space 3}        .
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode>cowgroupid             {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 41.17977{col 45}{space 2} 6.871314{col 73}{space 4} 29.69281{col 86}{space 3} 57.11059
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{txt}LR test vs. nbinomial model: {txt}chi2({res}2{txt}) ={res} 3401.72{col 59}{txt}Prob > chi2 ={res}{col 73}0.0000

{txt}{p 0 6 4 79}Note: {help j_mixedlr##|_new:LR test is conservative} and provided only for reference.{p_end}

{com}. est store mod1  
{txt}
{com}. coefplot (mod1, label(Model 1) pstyle(p5) msymbol(T)), drop(_cons) xline(0) xtitle(Average Marginal Effects) ///
>          order(v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///
>                Weighted_Conflict_Index_lag_log cw_binary log_rgdppc_lag muspct_lagged groupsize excl_groups_count) ///
>          coeflabel(v2exl_legitideolcr_0_lagged = "Leader Nationalism" xpolity2 = "Democracy" status_autonomy = "Regional Autonomy" status_powerless = "Powerless" status_discrim = "Discriminated" status_separatist = "Separatist" ///
>                            Weighted_Conflict_Index_lag_log = "Political Instability" cw_binary = "Civil War" log_rgdppc_lag = "Econ Development" muspct_lagged = "Muslims" groupsize = "Group Size" ///
>                                    excl_groups_count = "Terrorist Outbidding") ///
>          ciopts(recast(rcap))
{res}{txt}
{com}. estimates clear
{txt}
{com}. 
. 
. ***** Table 3 *****
. * Model 1, no mixed effects
. xtset cowgroupid year, yearly
{res}{txt}{col 8}panel variable:  {res}cowgroupid (unbalanced)
{txt}{col 9}time variable:  {res}{col 25}year, 1970 to 2009, but with gaps
{txt}{col 17}delta:  {res}1 year
{txt}
{com}. nbreg attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count ///                                                                                                           
>         if isrelevant==1 

{txt}Fitting Poisson model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-55243.525}  
Iteration 1:{space 3}log likelihood = {res:-54511.284}  
Iteration 2:{space 3}log likelihood = {res:-54501.185}  
Iteration 3:{space 3}log likelihood = {res:-54501.177}  
Iteration 4:{space 3}log likelihood = {res:-54501.177}  
{res}
{txt}Fitting constant-only model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-31878.726}  
Iteration 1:{space 3}log likelihood = {res:-8921.4575}  
Iteration 2:{space 3}log likelihood = {res:-8896.5713}  
Iteration 3:{space 3}log likelihood = {res:-8896.3489}  
Iteration 4:{space 3}log likelihood = {res:-8896.3489}  
{res}
{txt}Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-8700.2438}  (not concave)
Iteration 1:{space 3}log likelihood = {res:-8364.9658}  
Iteration 2:{space 3}log likelihood = {res:  -8174.96}  
Iteration 3:{space 3}log likelihood = {res: -8149.576}  
Iteration 4:{space 3}log likelihood = {res:-8149.2749}  
Iteration 5:{space 3}log likelihood = {res:-8149.2748}  
{res}
{txt}Negative binomial regression{col 49}Number of obs{col 67}= {res}    22,754
{txt}{col 49}LR chi2({res}12{txt}){col 67}= {res}   1494.15
{txt}{col 1}Dispersion{col 16}= {res}mean{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-8149.2748{txt}{col 49}Pseudo R2{col 67}= {res}    0.0840

{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2}  .430997{col 45}{space 2} .2043463{col 56}{space 1}    2.11{col 65}{space 3}0.035{col 73}{space 4} .0304855{col 86}{space 3} .8315084
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0875903{col 45}{space 2} .0138602{col 56}{space 1}    6.32{col 65}{space 3}0.000{col 73}{space 4} .0604248{col 86}{space 3} .1147557
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2}   2.2046{col 45}{space 2} .2073142{col 56}{space 1}   10.63{col 65}{space 3}0.000{col 73}{space 4} 1.798272{col 86}{space 3} 2.610929
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} 1.173185{col 45}{space 2} .1425995{col 56}{space 1}    8.23{col 65}{space 3}0.000{col 73}{space 4} .8936955{col 86}{space 3} 1.452675
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2} 1.650287{col 45}{space 2} .1483749{col 56}{space 1}   11.12{col 65}{space 3}0.000{col 73}{space 4} 1.359477{col 86}{space 3} 1.941096
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} 1.860951{col 45}{space 2}  .361435{col 56}{space 1}    5.15{col 65}{space 3}0.000{col 73}{space 4} 1.152551{col 86}{space 3}  2.56935
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2} .2451179{col 45}{space 2} .0134752{col 56}{space 1}   18.19{col 65}{space 3}0.000{col 73}{space 4} .2187069{col 86}{space 3} .2715289
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.412505{col 45}{space 2}  .107515{col 56}{space 1}   22.44{col 65}{space 3}0.000{col 73}{space 4}  2.20178{col 86}{space 3} 2.623231
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2}  .261451{col 45}{space 2} .0484673{col 56}{space 1}    5.39{col 65}{space 3}0.000{col 73}{space 4}  .166457{col 86}{space 3} .3564451
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2}-.0014603{col 45}{space 2} .0016576{col 56}{space 1}   -0.88{col 65}{space 3}0.378{col 73}{space 4}-.0047091{col 86}{space 3} .0017886
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.5866647{col 45}{space 2} .2416322{col 56}{space 1}   -2.43{col 65}{space 3}0.015{col 73}{space 4}-1.060255{col 86}{space 3}-.1130743
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0734406{col 45}{space 2} .0053543{col 56}{space 1}  -13.72{col 65}{space 3}0.000{col 73}{space 4}-.0839348{col 86}{space 3}-.0629464
{txt}{space 26}_cons {c |}{col 33}{res}{space 2} -6.33441{col 45}{space 2} .4135567{col 56}{space 1}  -15.32{col 65}{space 3}0.000{col 73}{space 4}-7.144967{col 86}{space 3}-5.523854
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
                       /lnalpha {c |}{col 33}{res}{space 2} 3.603834{col 45}{space 2} .0361324{col 73}{space 4} 3.533015{col 86}{space 3} 3.674652
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
                          alpha {c |}{col 33}{res}{space 2} 36.73881{col 45}{space 2} 1.327463{col 73}{space 4} 34.22702{col 86}{space 3} 39.43493
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
LR test of alpha=0: {help j_chibar##|_new:chibar2(01) = }{res}9.3e+04{col 56}{txt}Prob >= chibar2 = {res}0.000
{txt}
{com}. * Model 2, lagged DV
. menbreg attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count L.attack ///                                                                                                           
>         if isrelevant==1||cowcode: ||cowgroupid: , cov(uns) 
{res}{txt}
Fitting fixed-effects model:

Iteration 0:{space 3}log likelihood = {res:-15886.602}  
Iteration 1:{space 3}log likelihood = {res:-10650.722}  
Iteration 2:{space 3}log likelihood = {res:-8160.4858}  
Iteration 3:{space 3}log likelihood = {res:-7874.4895}  
Iteration 4:{space 3}log likelihood = {res:-7702.4141}  
Iteration 5:{space 3}log likelihood = {res:-7641.2951}  
Iteration 6:{space 3}log likelihood = {res:-7633.3976}  
Iteration 7:{space 3}log likelihood = {res:-7633.2922}  
Iteration 8:{space 3}log likelihood = {res:-7633.2921}  

Refining starting values:

Grid node 0:{space 3}log likelihood = {res:-7205.6807}

Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-7205.6807}  
Iteration 1:{space 3}log likelihood = {res:-6925.2458}  (backed up)
Iteration 2:{space 3}log likelihood = {res:-6512.3431}  
Iteration 3:{space 3}log likelihood = {res:-6492.5241}  
Iteration 4:{space 3}log likelihood = {res:-6474.5776}  
Iteration 5:{space 3}log likelihood = {res:-6466.1602}  
Iteration 6:{space 3}log likelihood = {res:-6458.1159}  
Iteration 7:{space 3}log likelihood = {res:-6450.4204}  
Iteration 8:{space 3}log likelihood = {res:-6436.1231}  
Iteration 9:{space 3}log likelihood = {res:-6337.9095}  
Iteration 10:{space 2}log likelihood = {res:-6293.5388}  
Iteration 11:{space 2}log likelihood = {res:-6278.7844}  
Iteration 12:{space 2}log likelihood = {res:-6275.0855}  
Iteration 13:{space 2}log likelihood = {res:-6272.7779}  
Iteration 14:{space 2}log likelihood = {res: -6271.489}  
Iteration 15:{space 2}log likelihood = {res:-6270.8703}  
Iteration 16:{space 2}log likelihood = {res:-6270.6315}  
Iteration 17:{space 2}log likelihood = {res: -6270.571}  
Iteration 18:{space 2}log likelihood = {res:-6270.5757}  
Iteration 19:{space 2}log likelihood = {res:-6270.5962}  
Iteration 20:{space 2}log likelihood = {res:-6270.6161}  
Iteration 21:{space 2}log likelihood = {res:-6270.6321}  
Iteration 22:{space 2}log likelihood = {res:-6270.6442}  
Iteration 23:{space 2}log likelihood = {res:-6270.6532}  
Iteration 24:{space 2}log likelihood = {res:-6270.6598}  
Iteration 25:{space 2}log likelihood = {res:-6270.6647}  
Iteration 26:{space 2}log likelihood = {res:-6270.6684}  
Iteration 27:{space 2}log likelihood = {res: -6270.671}  
Iteration 28:{space 2}log likelihood = {res:-6270.6729}  
Iteration 29:{space 2}log likelihood = {res:-6270.6743}  
Iteration 30:{space 2}log likelihood = {res:-6270.6754}  
Iteration 31:{space 2}log likelihood = {res:-6270.6761}  
Iteration 32:{space 2}log likelihood = {res:-6270.6767}  
Iteration 33:{space 2}log likelihood = {res:-6270.6771}  
{res}
{txt}Mixed-effects nbinomial regression{col 49}{txt}Number of obs{col 67}={res}{col 69}    22,033
{txt}Overdispersion: {col 28}{res}mean

{txt}{hline 16}{c TT}{hline 44}
{col 17}{txt}{c |}{col 23}No. of{col 36}Observations per Group
{col 2}{txt}Group Variable{col 17}{c |}{col 23}Groups{col 33}Minimum{col 44}Average{col 55}Maximum
{txt}{hline 16}{c +}{hline 44}
{col 9}{res}cowcode{col 17}{txt}{c |}{res}{col 21}     132{col 31}        4{col 42}    166.9{col 53}    1,803
{col 6}{res}cowgroupid{col 17}{txt}{c |}{res}{col 21}     756{col 31}        1{col 42}     29.1{col 53}       39
{txt}{hline 16}{c BT}{hline 44}

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}13{txt}){col 67}={res}{col 70}   632.93
{txt}Log likelihood = {res}-6270.6771{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2} 1.106752{col 45}{space 2} .2888838{col 56}{space 1}    3.83{col 65}{space 3}0.000{col 73}{space 4} .5405505{col 86}{space 3} 1.672954
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0283596{col 45}{space 2} .0165389{col 56}{space 1}    1.71{col 65}{space 3}0.086{col 73}{space 4}-.0040561{col 86}{space 3} .0607753
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2} .1468748{col 45}{space 2} .2815387{col 56}{space 1}    0.52{col 65}{space 3}0.602{col 73}{space 4}-.4049311{col 86}{space 3} .6986806
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2}-.0257073{col 45}{space 2} .2124099{col 56}{space 1}   -0.12{col 65}{space 3}0.904{col 73}{space 4} -.442023{col 86}{space 3} .3906084
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2}  .745896{col 45}{space 2} .2231265{col 56}{space 1}    3.34{col 65}{space 3}0.001{col 73}{space 4}  .308576{col 86}{space 3} 1.183216
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} .8851749{col 45}{space 2}  .402691{col 56}{space 1}    2.20{col 65}{space 3}0.028{col 73}{space 4}  .095915{col 86}{space 3} 1.674435
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2}  .082889{col 45}{space 2} .0132729{col 56}{space 1}    6.24{col 65}{space 3}0.000{col 73}{space 4} .0568746{col 86}{space 3} .1089033
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.169332{col 45}{space 2} .1254758{col 56}{space 1}   17.29{col 65}{space 3}0.000{col 73}{space 4} 1.923404{col 86}{space 3}  2.41526
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2}  .617918{col 45}{space 2} .1249465{col 56}{space 1}    4.95{col 65}{space 3}0.000{col 73}{space 4} .3730273{col 86}{space 3} .8628086
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2} .0105313{col 45}{space 2} .0053413{col 56}{space 1}    1.97{col 65}{space 3}0.049{col 73}{space 4} .0000626{col 86}{space 3} .0210001
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.7659899{col 45}{space 2} .7888616{col 56}{space 1}   -0.97{col 65}{space 3}0.332{col 73}{space 4} -2.31213{col 86}{space 3} .7801504
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0384549{col 45}{space 2} .0165888{col 56}{space 1}   -2.32{col 65}{space 3}0.020{col 73}{space 4}-.0709684{col 86}{space 3}-.0059413
{txt}{space 31} {c |}
{space 25}attack {c |}
{space 28}L1. {c |}{col 33}{res}{space 2} .0274223{col 45}{space 2} .0029518{col 56}{space 1}    9.29{col 65}{space 3}0.000{col 73}{space 4} .0216368{col 86}{space 3} .0332078
{txt}{space 31} {c |}
{space 26}_cons {c |}{col 33}{res}{space 2}-15.53859{col 45}{space 2}   1.3069{col 56}{space 1}  -11.89{col 65}{space 3}0.000{col 73}{space 4}-18.10007{col 86}{space 3}-12.97711
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 23}/lnalpha {c |}{col 33}{res}{space 2} 1.451193{col 45}{space 2} .0431569{col 56}{space 1}   33.63{col 65}{space 3}0.000{col 73}{space 4} 1.366607{col 86}{space 3} 1.535779
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode                        {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 3.83e-23{col 45}{space 2} 4.69e-12{col 73}{space 4}        .{col 86}{space 3}        .
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode>cowgroupid             {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 33.80028{col 45}{space 2} 5.771141{col 73}{space 4} 24.18711{col 86}{space 3} 47.23419
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{txt}LR test vs. nbinomial model: {txt}{help j_chibar##|_new:chibar2(01) =}{res} 2725.23{col 55}{txt}Prob >= chibar2 = {res}{col 73}0.0000
{txt}
{com}. * Model 3, alternative nationalism measure: nationalism_executive_1
. menbreg attack nationalism_executive_1 xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count ///                                                                                                           
>         if isrelevant==1||cowcode: ||cowgroupid: ,cov(uns)
{res}{txt}
Fitting fixed-effects model:

Iteration 0:{space 3}log likelihood = {res:-23915.125}  
Iteration 1:{space 3}log likelihood = {res:-12149.452}  
Iteration 2:{space 3}log likelihood = {res: -11625.49}  (not concave)
Iteration 3:{space 3}log likelihood = {res:-10246.722}  (not concave)
Iteration 4:{space 3}log likelihood = {res:-8237.1996}  (not concave)
Iteration 5:{space 3}log likelihood = {res:-8084.3554}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-7923.5431}  
Iteration 7:{space 3}log likelihood = {res:-7618.5801}  
Iteration 8:{space 3}log likelihood = {res:-7504.6718}  
Iteration 9:{space 3}log likelihood = {res:-7493.2906}  
Iteration 10:{space 2}log likelihood = {res:-7493.2173}  
Iteration 11:{space 2}log likelihood = {res:-7493.2173}  

Refining starting values:

Grid node 0:{space 3}log likelihood = {res:-7101.8338}

Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-7101.8338}  (not concave)
Iteration 1:{space 3}log likelihood = {res:-6169.6178}  
Iteration 2:{space 3}log likelihood = {res:-5975.4598}  
Iteration 3:{space 3}log likelihood = {res:-5903.8207}  
Iteration 4:{space 3}log likelihood = {res:-5877.2102}  
Iteration 5:{space 3}log likelihood = {res:-5869.7128}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-5869.2118}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-5869.7961}  (not concave)
Iteration 8:{space 3}log likelihood = {res:-5869.8163}  
Iteration 9:{space 3}log likelihood = {res:-5869.6844}  
Iteration 10:{space 2}log likelihood = {res:-5869.4488}  
Iteration 11:{space 2}log likelihood = {res:-5869.1906}  
Iteration 12:{space 2}log likelihood = {res: -5869.136}  
Iteration 13:{space 2}log likelihood = {res:-5869.1463}  
Iteration 14:{space 2}log likelihood = {res:-5869.1751}  
Iteration 15:{space 2}log likelihood = {res: -5869.202}  
Iteration 16:{space 2}log likelihood = {res: -5869.223}  
Iteration 17:{space 2}log likelihood = {res:-5869.2385}  
Iteration 18:{space 2}log likelihood = {res:-5869.2496}  
Iteration 19:{space 2}log likelihood = {res:-5869.2577}  
Iteration 20:{space 2}log likelihood = {res:-5869.2634}  
Iteration 21:{space 2}log likelihood = {res:-5869.2676}  
Iteration 22:{space 2}log likelihood = {res:-5869.2705}  
Iteration 23:{space 2}log likelihood = {res:-5869.2727}  
Iteration 24:{space 2}log likelihood = {res:-5869.2742}  
Iteration 25:{space 2}log likelihood = {res:-5869.2753}  
Iteration 26:{space 2}log likelihood = {res: -5869.276}  
Iteration 27:{space 2}log likelihood = {res:-5869.2766}  
Iteration 28:{space 2}log likelihood = {res: -5869.277}  
{res}
{txt}Mixed-effects nbinomial regression{col 49}{txt}Number of obs{col 67}={res}{col 69}    19,301
{txt}Overdispersion: {col 28}{res}mean

{txt}{hline 16}{c TT}{hline 44}
{col 17}{txt}{c |}{col 23}No. of{col 36}Observations per Group
{col 2}{txt}Group Variable{col 17}{c |}{col 23}Groups{col 33}Minimum{col 44}Average{col 55}Maximum
{txt}{hline 16}{c +}{hline 44}
{col 9}{res}cowcode{col 17}{txt}{c |}{res}{col 21}     127{col 31}        6{col 42}    152.0{col 53}    1,618
{col 6}{res}cowgroupid{col 17}{txt}{c |}{res}{col 21}     725{col 31}        1{col 42}     26.6{col 53}       34
{txt}{hline 16}{c BT}{hline 44}

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}12{txt}){col 67}={res}{col 70}   461.74
{txt}Log likelihood = {res} -5869.277{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}nationalism_executive_1 {c |}{col 33}{res}{space 2} .8289541{col 45}{space 2} .2694594{col 56}{space 1}    3.08{col 65}{space 3}0.002{col 73}{space 4} .3008233{col 86}{space 3} 1.357085
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0583627{col 45}{space 2} .0186895{col 56}{space 1}    3.12{col 65}{space 3}0.002{col 73}{space 4} .0217319{col 86}{space 3} .0949935
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2}  .509021{col 45}{space 2} .3062815{col 56}{space 1}    1.66{col 65}{space 3}0.097{col 73}{space 4}-.0912797{col 86}{space 3} 1.109322
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} .0664172{col 45}{space 2} .2197673{col 56}{space 1}    0.30{col 65}{space 3}0.762{col 73}{space 4}-.3643187{col 86}{space 3} .4971532
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2} 1.061575{col 45}{space 2} .2383105{col 56}{space 1}    4.45{col 65}{space 3}0.000{col 73}{space 4} .5944952{col 86}{space 3} 1.528655
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} 1.283814{col 45}{space 2} .4222146{col 56}{space 1}    3.04{col 65}{space 3}0.002{col 73}{space 4} .4562885{col 86}{space 3} 2.111339
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2} .0993587{col 45}{space 2} .0146246{col 56}{space 1}    6.79{col 65}{space 3}0.000{col 73}{space 4} .0706949{col 86}{space 3} .1280224
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.331799{col 45}{space 2}  .138009{col 56}{space 1}   16.90{col 65}{space 3}0.000{col 73}{space 4} 2.061306{col 86}{space 3} 2.602292
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2} .6028122{col 45}{space 2} .1336384{col 56}{space 1}    4.51{col 65}{space 3}0.000{col 73}{space 4} .3408859{col 86}{space 3} .8647386
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2}  .009765{col 45}{space 2} .0058624{col 56}{space 1}    1.67{col 65}{space 3}0.096{col 73}{space 4} -.001725{col 86}{space 3} .0212551
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.6667649{col 45}{space 2} .8112297{col 56}{space 1}   -0.82{col 65}{space 3}0.411{col 73}{space 4}-2.256746{col 86}{space 3} .9232161
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0498929{col 45}{space 2} .0161329{col 56}{space 1}   -3.09{col 65}{space 3}0.002{col 73}{space 4}-.0815129{col 86}{space 3}-.0182729
{txt}{space 26}_cons {c |}{col 33}{res}{space 2} -15.4355{col 45}{space 2} 1.355424{col 56}{space 1}  -11.39{col 65}{space 3}0.000{col 73}{space 4}-18.09208{col 86}{space 3}-12.77892
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 23}/lnalpha {c |}{col 33}{res}{space 2} 1.423659{col 45}{space 2} .0437851{col 56}{space 1}   32.51{col 65}{space 3}0.000{col 73}{space 4} 1.337842{col 86}{space 3} 1.509476
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode                        {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 4.24e-25{col 45}{space 2} 4.73e-13{col 73}{space 4}        .{col 86}{space 3}        .
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode>cowgroupid             {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 38.19133{col 45}{space 2} 6.345317{col 73}{space 4} 27.57665{col 86}{space 3} 52.89177
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{txt}LR test vs. nbinomial model: {txt}{help j_chibar##|_new:chibar2(01) =}{res} 3247.88{col 55}{txt}Prob >= chibar2 = {res}{col 73}0.0000
{txt}
{com}. * Model 4, region
. sort cowcode year
{txt}
{com}. replace americas = americas[_n-1] if missing(americas) 
{txt}(1,024 real changes made)

{com}. replace europe = europe[_n-1] if missing(europe) 
{txt}(1,024 real changes made)

{com}. replace africa = africa[_n-1] if missing(africa) 
{txt}(1,024 real changes made)

{com}. replace middleeast = middleeast[_n-1] if missing(middleeast) 
{txt}(1,024 real changes made)

{com}. replace asia = asia[_n-1] if missing(asia) 
{txt}(0 real changes made)

{com}. replace oceanania = oceanania[_n-1] if missing(oceanania) 
{txt}(1,024 real changes made)

{com}. menbreg attack v2exl_legitideolcr_0_lagged xpolity2 status_autonomy status_powerless status_discrim status_separatist ///        
>     Weighted_Conflict_Index_lag_log cw_binary ///                                                                   
>         log_rgdppc_lag ///                                                                                                   
>         muspct_lagged ///                                                                                                                         
>         groupsize excl_groups_count americas europe africa middleeast asia ///                                                                                                           
>         if isrelevant==1||cowcode: ||cowgroupid: ,cov(uns)
{res}{txt}
Fitting fixed-effects model:

Iteration 0:{space 3}log likelihood = {res:-26657.347}  
Iteration 1:{space 3}log likelihood = {res:-15183.991}  (not concave)
Iteration 2:{space 3}log likelihood = {res:-10663.097}  (not concave)
Iteration 3:{space 3}log likelihood = {res:-9364.7938}  (not concave)
Iteration 4:{space 3}log likelihood = {res:  -9119.11}  (not concave)
Iteration 5:{space 3}log likelihood = {res:-8777.9065}  (not concave)
Iteration 6:{space 3}log likelihood = {res:-8586.0622}  (not concave)
Iteration 7:{space 3}log likelihood = {res:-8402.2408}  
Iteration 8:{space 3}log likelihood = {res:-8182.1211}  
Iteration 9:{space 3}log likelihood = {res:-8177.5633}  
Iteration 10:{space 2}log likelihood = {res:-8120.3975}  
Iteration 11:{space 2}log likelihood = {res:-8116.1859}  
Iteration 12:{space 2}log likelihood = {res:-8116.1446}  
Iteration 13:{space 2}log likelihood = {res:-8116.1445}  

Refining starting values:

Grid node 0:{space 3}log likelihood = {res:-7674.6121}

Fitting full model:
{res}
{txt}Iteration 0:{space 3}log likelihood = {res:-7674.6121}  
Iteration 1:{space 3}log likelihood = {res:-6694.1091}  (backed up)
Iteration 2:{space 3}log likelihood = {res:-6537.8083}  
Iteration 3:{space 3}log likelihood = {res:-6472.5572}  
Iteration 4:{space 3}log likelihood = {res:-6450.1458}  
Iteration 5:{space 3}log likelihood = {res:-6448.3747}  
Iteration 6:{space 3}log likelihood = {res:-6446.7952}  
Iteration 7:{space 3}log likelihood = {res:-6445.3916}  
Iteration 8:{space 3}log likelihood = {res:-6444.3218}  
Iteration 9:{space 3}log likelihood = {res:-6442.9858}  
Iteration 10:{space 2}log likelihood = {res:-6442.3892}  
Iteration 11:{space 2}log likelihood = {res:-6442.2072}  
Iteration 12:{space 2}log likelihood = {res:-6442.2146}  
Iteration 13:{space 2}log likelihood = {res:-6442.2842}  
Iteration 14:{space 2}log likelihood = {res:-6442.3603}  
Iteration 15:{space 2}log likelihood = {res:-6442.4258}  
Iteration 16:{space 2}log likelihood = {res:-6442.4784}  
Iteration 17:{space 2}log likelihood = {res:-6442.5199}  
Iteration 18:{space 2}log likelihood = {res:-6442.5522}  
Iteration 19:{space 2}log likelihood = {res:-6442.5646}  
Iteration 20:{space 2}log likelihood = {res:-6442.5757}  
Iteration 21:{space 2}log likelihood = {res:-6442.5855}  
Iteration 22:{space 2}log likelihood = {res:-6442.5943}  
Iteration 23:{space 2}log likelihood = {res:-6442.6022}  
Iteration 24:{space 2}log likelihood = {res:-6442.6089}  
Iteration 25:{space 2}log likelihood = {res:-6442.6151}  
Iteration 26:{space 2}log likelihood = {res:-6442.6204}  
Iteration 27:{space 2}log likelihood = {res:-6442.6251}  
Iteration 28:{space 2}log likelihood = {res:-6442.6293}  
Iteration 29:{space 2}log likelihood = {res:-6442.6329}  
Iteration 30:{space 2}log likelihood = {res:-6442.6395}  
Iteration 31:{space 2}log likelihood = {res:-6442.6419}  
Iteration 32:{space 2}log likelihood = {res:-6442.6507}  
Iteration 33:{space 2}log likelihood = {res:-6442.6511}  
Iteration 34:{space 2}log likelihood = {res:-6442.6554}  
Iteration 35:{space 2}log likelihood = {res:-6442.6575}  
Iteration 36:{space 2}log likelihood = {res:-6442.6586}  
Iteration 37:{space 2}log likelihood = {res:-6442.6586}  (backed up)
{res}
{txt}Mixed-effects nbinomial regression{col 49}{txt}Number of obs{col 67}={res}{col 69}    22,754
{txt}Overdispersion: {col 28}{res}mean

{txt}{hline 16}{c TT}{hline 44}
{col 17}{txt}{c |}{col 23}No. of{col 36}Observations per Group
{col 2}{txt}Group Variable{col 17}{c |}{col 23}Groups{col 33}Minimum{col 44}Average{col 55}Maximum
{txt}{hline 16}{c +}{hline 44}
{col 9}{res}cowcode{col 17}{txt}{c |}{res}{col 21}     132{col 31}        6{col 42}    172.4{col 53}    1,899
{col 6}{res}cowgroupid{col 17}{txt}{c |}{res}{col 21}     757{col 31}        1{col 42}     30.1{col 53}       40
{txt}{hline 16}{c BT}{hline 44}

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}17{txt}){col 67}={res}{col 70}   541.38
{txt}Log likelihood = {res}-6442.6586{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{hline 32}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}                         attack{col 33}{c |}      Coef.{col 45}   Std. Err.{col 57}      z{col 65}   P>|z|{col 73}     [95% Con{col 86}f. Interval]
{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}v2exl_legitideolcr_0_lagged {c |}{col 33}{res}{space 2} 1.028264{col 45}{space 2} .3131316{col 56}{space 1}    3.28{col 65}{space 3}0.001{col 73}{space 4} .4145369{col 86}{space 3}  1.64199
{txt}{space 23}xpolity2 {c |}{col 33}{res}{space 2} .0472135{col 45}{space 2} .0180176{col 56}{space 1}    2.62{col 65}{space 3}0.009{col 73}{space 4} .0118997{col 86}{space 3} .0825273
{txt}{space 16}status_autonomy {c |}{col 33}{res}{space 2} .3866999{col 45}{space 2} .3042592{col 56}{space 1}    1.27{col 65}{space 3}0.204{col 73}{space 4}-.2096372{col 86}{space 3} .9830369
{txt}{space 15}status_powerless {c |}{col 33}{res}{space 2} .1259065{col 45}{space 2} .2234128{col 56}{space 1}    0.56{col 65}{space 3}0.573{col 73}{space 4}-.3119744{col 86}{space 3} .5637875
{txt}{space 17}status_discrim {c |}{col 33}{res}{space 2}  1.04727{col 45}{space 2} .2376309{col 56}{space 1}    4.41{col 65}{space 3}0.000{col 73}{space 4} .5815225{col 86}{space 3} 1.513018
{txt}{space 14}status_separatist {c |}{col 33}{res}{space 2} 1.265281{col 45}{space 2} .4344753{col 56}{space 1}    2.91{col 65}{space 3}0.004{col 73}{space 4} .4137247{col 86}{space 3} 2.116837
{txt}Weighted_Conflict_Index_lag_log {c |}{col 33}{res}{space 2}  .109701{col 45}{space 2} .0142186{col 56}{space 1}    7.72{col 65}{space 3}0.000{col 73}{space 4} .0818332{col 86}{space 3} .1375689
{txt}{space 22}cw_binary {c |}{col 33}{res}{space 2} 2.445385{col 45}{space 2} .1342016{col 56}{space 1}   18.22{col 65}{space 3}0.000{col 73}{space 4} 2.182355{col 86}{space 3} 2.708416
{txt}{space 17}log_rgdppc_lag {c |}{col 33}{res}{space 2} .9369067{col 45}{space 2} .1563319{col 56}{space 1}    5.99{col 65}{space 3}0.000{col 73}{space 4} .6305018{col 86}{space 3} 1.243312
{txt}{space 18}muspct_lagged {c |}{col 33}{res}{space 2} .0072644{col 45}{space 2} .0070071{col 56}{space 1}    1.04{col 65}{space 3}0.300{col 73}{space 4}-.0064693{col 86}{space 3} .0209981
{txt}{space 22}groupsize {c |}{col 33}{res}{space 2}-.0964199{col 45}{space 2} .8549401{col 56}{space 1}   -0.11{col 65}{space 3}0.910{col 73}{space 4}-1.772072{col 86}{space 3} 1.579232
{txt}{space 14}excl_groups_count {c |}{col 33}{res}{space 2}-.0414634{col 45}{space 2} .0177064{col 56}{space 1}   -2.34{col 65}{space 3}0.019{col 73}{space 4}-.0761674{col 86}{space 3}-.0067594
{txt}{space 23}americas {c |}{col 33}{res}{space 2} .8819805{col 45}{space 2}  2.39772{col 56}{space 1}    0.37{col 65}{space 3}0.713{col 73}{space 4}-3.817465{col 86}{space 3} 5.581426
{txt}{space 25}europe {c |}{col 33}{res}{space 2} 1.437795{col 45}{space 2} 2.319321{col 56}{space 1}    0.62{col 65}{space 3}0.535{col 73}{space 4}-3.107991{col 86}{space 3}  5.98358
{txt}{space 25}africa {c |}{col 33}{res}{space 2} 3.143876{col 45}{space 2} 2.333032{col 56}{space 1}    1.35{col 65}{space 3}0.178{col 73}{space 4}-1.428783{col 86}{space 3} 7.716535
{txt}{space 21}middleeast {c |}{col 33}{res}{space 2} 2.515473{col 45}{space 2} 2.392666{col 56}{space 1}    1.05{col 65}{space 3}0.293{col 73}{space 4}-2.174067{col 86}{space 3} 7.205013
{txt}{space 27}asia {c |}{col 33}{res}{space 2} 3.201054{col 45}{space 2} 2.334286{col 56}{space 1}    1.37{col 65}{space 3}0.170{col 73}{space 4}-1.374064{col 86}{space 3} 7.776171
{txt}{space 26}_cons {c |}{col 33}{res}{space 2}-21.36952{col 45}{space 2} 2.786072{col 56}{space 1}   -7.67{col 65}{space 3}0.000{col 73}{space 4}-26.83012{col 86}{space 3}-15.90892
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 23}/lnalpha {c |}{col 33}{res}{space 2} 1.604745{col 45}{space 2} .0415452{col 56}{space 1}   38.63{col 65}{space 3}0.000{col 73}{space 4} 1.523318{col 86}{space 3} 1.686172
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode                        {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 3.18e-28{col 45}{space 2} 2.28e-14{col 73}{space 4}        .{col 86}{space 3}        .
{txt}{hline 32}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}cowcode>cowgroupid             {col 33}{c |}
                      var(_cons){c |}{col 33}{res}{space 2} 39.37776{col 45}{space 2} 6.610732{col 73}{space 4} 28.33684{col 86}{space 3} 54.72056
{txt}{hline 32}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{txt}LR test vs. nbinomial model: {txt}{help j_chibar##|_new:chibar2(01) =}{res} 3346.97{col 55}{txt}Prob >= chibar2 = {res}{col 73}0.0000
{txt}
{com}. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
. 
{txt}end of do-file

{com}. exit, clear
